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Ph; abstract 

Q . We present deep Ka-band {v k, 33 GHz) observations of the CMB made with the 

$-H ' extended Very Small Array (VSA). This configuration produces a naturally weighted 

\ synthesized FWHM beamwidth of ~ 11 arcmin which covers an £-range of 300 to 1500. 

■ On these scales, foreground extragalactic sources can be a major source of contamina- 

^ ' tion to the CMB anisotropy. This problem has been alleviated by identifying sources 

• i-H , at 15 GHz with the Ryle Telescope and then monitoring these sources at 33 GHz using 

' a single baseline interferometer co-located with the VSA. Sources with flux densities 

^20 mJy at 33 GHz are subtracted from the data. In addition, we calculate a statis- 
tical correction for the small residual contribution from weaker sources that are below 
the detection limit of the survey. 

The CMB power spectrum corrected for Galactic foregrounds and extragalactic 
point sources is presented. A total £-range of 150 — 1500 is achieved by combining the 
complete extended array data with earlier VSA data in a compact configuration. Our 
resolution of At k, 60 allows the first 3 acoustic peaks to be clearly delineated. The 
is achieved by using mosaiced observations in 7 regions covering a total area of 82 sq. 
degrees. There is good agreement with WMAP data up to £ = 700 where WMAP data 
run out of resolution. For higher ^-values out to £ = 1500, the agreement in power 
spectrum amplitudes with other experiments is also very good despite differences in 
frequency and observing technique. 

Key words: cosmology:observations - cosmic microwave background - techniques: 
interferometric 
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I INTRODUCTION 

The angular power spectrum of primordial anisotropies in 
the cosmic microwave background (CMB) has become an 

II Present address: Instituto de Astrofisica de Canarias, 38200 La 
Laguna, Tenerife, Spain 
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important tool in the era of "precision cosmology". Since 
the first statistical detection of CMB fluctuations on large 
angular scales = 2 - 30) by the COBE-'DMR instru- 
ment (Smoot et al. 1992), several experiments have de- 
tected acoustic peaks in the power spectrum in the ^-range 
100-1000 (Lee et al. 2001; Netterfield et al. 2002; Halverson 
et al. 2002; Benoit et al. 2003; Scott et al. 2003) and a fall-off 
in power at high ^-values (Dawson et al. 2002; Grainge et 
al. 2003; Kuo et al. 2004; Readhead et al. 2004). More re- 
cently, the Wilkinson Microwave Anisotropy Probe, hence- 
forth WMAP, has provided unprecedented measurements 
over the l-range 2 — 700 (Bennett et al. 2003a; Hinshaw et al. 
2003a). The WMAP 1-year power spectrum is cosmic vari- 
ance limited up to ^ = 350 and delineates the first 2 peaks at 
£ ~ 220 and £ ~ 550 with exceptional signal-to-noise ratios. 
The new data have provided detailed cosmological informa- 
tion on a wide range of parameters (Spergel et al. 2003) and 
have raised new questions to be answered. Ifowever, the an- 
gular resolution of WMAP limits the power spectrum to 
£ ^ 800. At high ^-values, the origin of the excess power ini- 
tially observed by the Cosmic Background Imager (CBI) at 
£ = 2000-4000 (Mason et al. 2003) has also generated much 
interest. It is clear that the high-^ CMB power spectrum is 
one of the challenges for future CMB experiments including 
the ESA Planck satellite (Tauber 2001; Lawrence 2003) due 
for launch in 2007. 

The Very Small Array, henceforth the VSA (Watson et 
al. 2003, hereafter Paper I), is a purpose-built radio inter- 
ferometer that has measured the CMB angular power spec- 
trum between £ = 150 and 900 in a compact array config- 
uration (Scott et al. 2003, Paper III) and more recently up 
to ^ = 1400 in an extended array (Grainge et al. 2003). This 
paper describes the complete set of extended array data ob- 
served during October 2001 — June 2003. The data cover 
a larger area of sky than those analysed by Grainge et al. 
(2003) by adding more pointings to the previous VSA fields 
and by observing 4 new fields. We use mosaicing techniques 
to increase the sensitivity, reduce the sample variance and 
to facilitate finer £ resolution — or equivalently reduce cor- 
relations between £ bins — over the range £ — 300 — 1500. 

The paper is organised as follows. Section|5|summarises 
the telescope parameters and the extended array configu- 
ration. In section 13 we describe the observations, data re- 
duction and calibration of the data, including a range of 
data checks. In section^Jwe discuss the various foregrounds, 
particularly discrete radio sources, and the corrections that 
were made to the data. The main results, CMB mosaic maps 
and CMB power spectrum covering £ — 150 — 1500, are pre- 
sented in section |S] A morphological analysis of the power 
spectrum is given in section 15.31 Section |S| is a discussion 
of the results and comparisons with other data, followed by 
conclusions in section |7| The cosmological interpretation of 
these data, both on their own and combined with other data, 
is described in a companion paper (Rebolo et al. 2004). 



2 THE VSA EXTENDED ARRAY 

The VSA is a 14-element interferometer operating in the Ka- 
band (26 — 36 GHz) situated at the high and dry site of the 
El Teide observatory, in Tenerife, at an altitude of 2340 m 
(see Paper I). Each antenna comprises a conical corrugated 



horn feeding a paraboloid mirror and mounted anywhere 
on a 4-mx3-m tip-tilt table located in a metal enclosure to 
minimise ground emission. Due to the geometry of the ta- 
ble and enclosure, the VSA declination range is restricted to 
—5° to -1-60°. In the "compact" configuration (Taylor et al. 
2003, hereafter Paper II), the mirrors were 143- mm in diam- 
eter giving a primary beam of 4°.6 FWHM at 34 GHz. The 
"extended array" , described here, has larger 322-mm diam- 
eter apertures allowing longer baselines while maintaining 
a high filling factor and therefore high temperature sensi- 
tivity. The longer baselines increase the angular resolution 
while the larger horns give a primary beam of 2°.l FWHM at 
33 GHz. This configuration has a total of 91 baselines with 
lengths ranging from 0.6 m to 2.5 m as shown in Fig.0 The 
maximum possible baseline length, set by the extent of the 
main tip-tilt table, is considerably higher with lengths up 
to ~ 4 m. However, for the data presented here, we chose 
a limited range of lengths to maximise the overall tempera- 
ture sensitivity of the array^ . This still allows a wide range 
of angular scales to be measured but with well-sampled u, v 
coverage (Fig. 0, which in turn provides good sampling of 
the power spectrum in £ and mapping with high fidelity. 
The filling factor of this configuration is ~ 1.6 times greater 
than that of the compact array, leading to an increase in the 
overall temperature sensitivity of the array. 

The VSA has so far operated in a single channel with 
instantaneous band-width of 1.5 GHz. We chose to use the 
higher end of the band (~ 33 GHz) to minimise foregrounds 
since the minimum contamination for total-power measure- 
ments is at ~ 70 GHz (Banday et al. 2003; Bennett et al. 
2003b). With an average system temperature of ~ 35 K, the 
VSA achieves an overall instantaneous point source sensitiv- 
ity of ~ 6 Jy s^^''^. This corresponds to a temperature sen- 
sitivity, over a synthesized beam area (ilsynth ~ 1 x 10"^ sr) 
of ~ 15 mK s"^''^. Note that the exact conversion from flux 
density to temperature depends on the beam area and on the 
u, V coverage which in turns depends on the declination and 
the flagging/filtering of the visibility data. A typical VSA 
field at declination 40° gives a total ^-range of ~ 300 — 1500 
and a naturally weighted synthesized beam of FWHM ~ 11 
arcmin over a ~ 2°.l field-of-view. The specifications for the 
VSA extended array are summarised in Table 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Observations 

The observations presented in this paper were made during 
the period October 2001- July 2003 with the array configu- 
ration described in section |5| They consist of a total of 33 
pointings that make up three 7-field mosaics (VSAl, VSA2, 
VSA3) and four 3-field mosaics (VSA5,VSA6,VSA7,VSA8) 
giving a total area coverage of 82 sq. degrees. The sky po- 
sitions for these 7 regions are depicted in Fig. |5| superposed 
on the 100 jim all-sky map based on DIRBE/IRAS data 
(Schlegel, Finkbeiner & Davis (1998), henceforth SFD98) 
in Galactic coordinates. The field centres and approximate 

^ The maximum baseline lengths possible with the current tip- 
tilt table, corresponding to £ ^ 2500, will be utilised in future 
CMB observations with the VSA. 
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Figure 1. Left: The extended array configuration of 14 antennas 
on tile tip-tiit table. Right: Tlie corresponding u, v coverage (units 
of wavelengtlis) for a 5 liour observation at declination +40° . 



Table 1. Specifications of the VSA in the extended array config- 
uration. 



Location 
Latitude 

Declination range 

No. of antennas (baselines) 

Range of baseline lengths 

Centre frequencies 

Bandwidth, Au 

System temperature, Tsys (K) 

Mirror diameters 

Primary beam 

Synthesized beam (FWHM) 

Range of £ 

Point source flux sensitivity 
Temperature sensitivity 



Izana, Tenerife (2340 m) 
+28° 18™ 
-5° <Doc < +60° 
14 (91) 
0.6 m —2.5 m 
33, 34 GHz 
1.5 GHz 
~ 35 K 
322 mm 
2°.l FWHM at 33 GHz 
11 arcmin 
- 300 - 1500 
~6 Jy s-i/2 

~ 15 mK s-i/2 









It 



300 SCO 100 

Galactic longitude (deg.) 

Figure 2. All-sky 100 ^tm map (SFD98) in Galactic coordinates. 
The grey-scale is linear from to 20 MJy sr~^. The locations of 
the 7 VSA regions are shown as circles, stretched in the horizontal 
direction due to the map projection, with diameters scaled up for 
clarity. The larger ellipses are the 7-field regions, labelled 1,2,3 
and the smaller ellipses, labelled 5,6,7,8, are the 3-field regions. 
The dashed and dotted lines indicate the declination limits for 
the VSA at +60° and -5° respectively. 



The results for the first 3 of the 7 fields in the VSAl, 
VSA2 and VSA3 regions (fields E,F and G) were pubhshed 
earlier (Grainge et al. 2003; Slosar et al. 2003). In this pa- 
per, we extend these regions (VSAl, VSA2 and VSA3) with 
a further 4 pointings per region (H,J,K and L) plus a further 
4 new regions (VSA5, VSA6, VSA7 and VSA8) each with 3 
pointings (E,F and G). This corresponds to a factor of ~ 4 in 
the amount of extended array data and a significant increase 
in sensitivity over the previous results. Furthermore, the in- 
crease in sky-coverage and mosaicing in each field allows an 
improvement in ^-resolution, or equivalently reduced bin-bin 
correlations, using mosaicing techniques (section 15.21 . 



integration times are listed in Table |5] The total effective 
integration time is « 6000 hours (250 days) after filtering 
and flagging of the data. 

The fields were carefully chosen to minimise Galactic 
and extragalactic emission. This includes avoiding bright 
galaxy clusters based on existing catalogues (Ebeling et al. 
1998; Abell 1958) and bright radio sources (> 500 mjy) 
based on the NVSS 1.4 GHz survey (Condon et al. 1998) 
and GB6 survey at 4.85 GHz (Gregory et al. 1996). Avoid- 
ing galaxy clusters is important due to the potential for 
Sunyaev-Zeldovich Effect (SZE) decrements to contaminate 
VSA data (see section !!^ . Diffuse Galactic emission is min- 
imised by choosing fields with low emission as predicted by 
templates of synchrotron, free-free and dust emission (see 
section [4.21 : all fields are at Galactic latitudes |&| It 27°. 
The VSA regions are distributed evenly in Right Ascension 
to optimise for 24 hour observations. Each field is observed 
on a daily basis for up to a maximum of 5 hours. Where pos- 
sible, observations are centred on transit to minimise ground 
spillover and system temperatures, giving a maximum Hour 
Angle range of ±2.5 hours. Calibrator observations of sev- 
eral bright radio sources were interleaved between fields for 
typically 10 min — 1 hour to provide fiux and phase calibra- 
tions (section 13.311 . At the same time as the VSA, a single- 
baseline interferometer monitors the brightest sources in the 
same patch of sky as part of the source-subtraction process 
described in more detail in section im 



3.2 Data reduction pipeline 

The data reduction and calibration procedures are similar 
to those used for the compact array data (Paper H). Each 
observation is analysed individually using the reduce pack- 
age written specifically for the VSA. The data reduction 
procedure is now highly developed; much of the correcting, 
flagging, filtering and re- weighting of the data are performed 
automatically. This results in data that are often close to the 
quality and sensitivity required for power spectrum estima- 
tion. We did however, identify residual atmospheric contam- 
ination which was not fully accounted for by the standard 
routines. Such low-level emission can easily be over-looked 
once stacked into large data sets with mostly good quality 
data. Fortunately, interferometers are much less sensitive to 
such signals than total-power instruments, since they "re- 
solve out" large angular scale emission due to incomplete 
u, V coverage at low spatial frequencies. The automatic rou- 
tines dealt with correlated emission either by discarding or 
down- weighting of noisy or contaminated data. Nevertheless, 
in a few rare instances, some low-level contamination still 
remained. They were initially detected by non-Gaussianity 
tests (section l3.5t which are very sensitive to such signals. 
We found that the best way to identify low-level contam- 
inated data was to examine the raw data by eye with no 
smoothing or calibration applied. Once they were identifled, 
we could remove all residual atmospheric emission or other 



4 C. Dickinson et al. 



Table 2. Celestial (J2000) coordinates for each of the 33 fields ob- 
served by the VSA in the extended array. The effective integration 
time is calculated after flagging and filtering of the data. VSAl, 
VSA2 and VSA3 are 7-field mosaics and VSA5, VSA6, VSA7 and 
VSA8 are 3-field mosaics giving a total of 33 pointings. 
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non-astronomical signals to well below the noise on time- 
scales of a few seconds. The contribution therefore to the 
overall data set is negligible. We also searched for noise cor- 
relations introduced by such emission. We show that such 
bias is negligible in section 1531 

We continue to use Fourier filtering to remove the ma- 
jority of local correlated signals including the so-called "spu- 
rious signal" described in Paper I. The longer baselines and 
larger horns of the extended array have significantly reduced 
this problem. The filtering removes typically 10 — 20 per cent 
of the data. The same fringe-rate filtering technique is also 
applied to the Sun and Moon. Data are filtered if the Sun 
and Moon are within 27° and 18° respectively, while if the 
Sun or Moon are within 9° of the field centre then the entire 
observation is fiagged. No residual Sun or Moon contamina- 
tion was detected after stacking the data typically integrated 
over 50 — 100 days. The data are then further smoothed by a 
factor of 4 to give 64 sec samples and a correction (discussed 



in section is applied for the atmospheric opacity which 
is typically a few per cent at hour angles of ~ 2.5 hours. 

The final step for each observation is the re-weighting 
of the data based on the r.m.s. noise of each baseline. This is 
important since the noise figures can vary substantially be- 
tween baselines. The re-weighting also allows the optimum 
overall noise level to be achieved. The data are then stacked 
together either in hour angle, or in the u, v plane, to reduce 
the noise level and to facilitate further flagging of any re- 
maining low-level bad data. Various consistency checks are 
performed, to verify the quality of data reduction; these are 
discussed in section 1531 Each VSA field contains ~ 10^ vis- 
ibilities, each of 64 sec integration, which are used directly 
to make maps (section l5.H . For power spectrum estimation 
(section l5.2ll . the data are binned in the u, v plane to reduce 
the number of data points. Each visibility has an associated 
weight calculated by the reduction pipeline by accounting 
for integration times, flagging and noise. The reliability of 
these weights is verified by comparing them with those cal- 
culated from the data themselves based on the scatter in 
each u, v cell containing several independent 64 sec sam- 
ples. Both approaches gave results consistent to within a 
few per cent. The long-term noise characteristics of the data 
are agree very well with the short-term estimates from indi- 
vidual visibilities. 

3.3 Amplitude and phase calibration 

Calibration of VSA data can be split into 4 parts; geomet- 
ric, amplitude, phase and atmospheric calibrations. The first 
stage is to calculate the exact geometry of the array. This 
requires knowledge of the x,y,z positions for each horn ac- 
curate to < 0.1 mm. This also requires calculating the other 
telescope parameters — including amplitudes, phases and 
the mean IF frequency — to converge on the correct ge- 
ometry, a total of ~ 400 parameters. We use a maximum- 
likelihood method developed specifically for interferometric 
arrays, introduced by Maisinger et al. (2003). The geometry 
is re-calculated from several long observations of Tau-A each 
time the array is modified. Thorough testing of the geome- 
try is made by calibrating observations of other bright radio 
sources. 

Amplitude and phase corrections are calculated from a 
single calibrator for each of the 91 complex channels (base- 
lines). An unresolved, non- variable bright radio source al- 
lows the measured fringes to be corrected for amplitude and 
phase. Table|5|hsts the primary calibration sources used by 
the VSA with their assumed fiux density 5* (Jy), or bright- 
ness temperature (K), and fiux density spectral index, a. 
Note that we adopt the convention S oc for the fiux 
density spectral index. The absolute calibration from cor- 
relator units to fiux density units (Jy) is described in sec- 
tion |5|1] The increased fiux sensitivity of the extended ar- 
ray, compared to the compact array, means that calibra- 
tions are reliable even for short (< 1 hour) observations. 
By cross-calibrating various bright radio sources with other 
sources, with no model-fitting, self-calibration or other post- 
processing techniques, we found that the majority of VSA 
visibility data have phases errors less than ~ 10°. We tested 
"secondary" (antenna-based) corrections based on calibra- 
tions made at the beginning and end of each observation. If 
necessary, this can correct for phase drifts during long ob- 
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Table 3 . Calibrator sources used for amplitude and phase calibra- 
tion. The absolute temperature scale is tied to a Jupiter bright- 
ness temperature of (146.6 ± 2.0) K at 33.0 GHz. Flux density 
spectral indices are in the convention S oc u~'^. 



Source RA/Dec Flux density (Jy) / Spectral 

(J2000) temperature (K) index a 

Tau-A 05'^35™ -f 22°01™ 333 Jy (34 GHz) -f0.30 

Cas-A 23^23™ + 58°49™ 168 Jy (34 GHz) -1-0.83 

Jupiter - 146.6 K (33 GHz) -2.24 



is consistent with the Mason et al. (1999) value at the la 
level. However, it corresponds to a considerable 8 per cent 
reduction in the VSA power spectrum values compared to 
previous VSA data. With that in mind, we also calculated 
the flux density of Tau-A and Cas-A from the 1-year WMAP 
data. The flux density values obtained from WMAP alone 
were found to be « 4 per cent lower than those using the 
original VSA value in good agreement with the new cali- 
bration. Earlier VSA power spectrum measurements should 
therefore be scaled by a factor of 0.92 to be consistent with 
the new calibration scheme used here. A recent paper de- 
scribing new CMB results from the CBI (Readhead et al. 
2004) derived a similar scaling factor of 0.94 for their data, 
based on the Jupiter temperature from WMAP. 

With angular extents of ~ 5 arcmin, both Tau-A and 
Cas-A are partially resolved by up to ~ 15 per cent on the 
longest baselines in the extended array. To correct for this, 
we used higher resolution maps, made at different wave- 
lengths, as a basis of a model. For Cas-A we used 32 GHz 
data taken with the 100 m Effelsberg telescope (W. Reich, 
priv. comm.). For Tau-A, a 1.4 GHz map made with the 
VLA (Bietenholz, Frail & Hester, 2001) was used and veri- 
fied with 850 /im map from the SCUBA instrument (Green 
2002). We adopt the flux density spectral indices given from 
Mason et al. (1999) (Table |^ to extrapolate between the 
frequencies. In a few rare instances, we also used Cyg-A 
(33.2 Jy at 34 GHz), Saturn (T = 135.4 K at 32 GHz) or 
Venus (T = 442.6 K at 32 GHz). 

3.5 Data Checks 

The VSA data are checked in a number of ways to ensure 
that the data are reduced correctly and that there are no 
systematics or contamination of the data. The usual checks 
consist of (i) inspecting the data by eye, (ii) parallel and in- 
dependent data reduction, (iii) stacking the data in different 
ways, (iv) data splits, (v) simulations and (vi) a variety of 
tests to search for residual non-Gaussianity. 

The data from each observation are displayed and ex- 
amined for contaminating signals such as correlated emission 
from atmospheric water clouds unresolved by the telescope 
beam. The worst periods of bad weather are easily identi- 
fied by an increase in system temperature and overall noise 
figures. However, unresolved clouds can result in correlated 
visibilities without a noticeable increase in system tempera- 
ture. We found that in a small portion of the data, correlated 
emission remained (section 13.211 . Although such data made 
little difference to the final noise level, there is a concern of 
introducing noise correlations which could affect the power 
spectrum in a biased fashion. Both the reduction software 
and the power spectrum estimate software (section 15.21 can 
deal with the full noise covariance matrix, but we assume 
that the off-diagonal elements were negligible to speed up 
the calculation of the power spectrum. We tested this as- 
sumption by calculating the non-diagonal elements of the 
full noise covariance matrix. On typical days the raw data 
are consistent with having zero correlations in the noise. 
We found that on the worst affected data, the automatic 
routines (weighting,filtering, atmospheric calibration and re- 
weighting) reduced the correlations to a very large degree. 
The off-diagonal elements were found to be negligible in vir- 
tually all cases. It is worth noting, that even if non-zero 



servations. However, in general, no correction was required 
since the VSA typically remained phase stable for periods 
of several days as demonstrated in Paper I. Uncorrelated 
phase errors of ~ 10° on individual visibilities have a negli- 
gible effect on CMB data since they tend to average to zero 
when combining a large number of visibilities. Also, the rel- 
atively low signal-to-noise ratios in the visibility data means 
that moderate phase errors (~ 10° or more) can be tolerated 
since we require a relatively low dynamic range. 

Data recorded before 17 July 2002 were taken at a 
frequency of 34 GHz. Due to a technical problem, we re- 
placed the tunable local oscillator (LO) with a fixed fre- 
quency quartz LO. This resulted in a new centre frequency 
of 33.0 GHz for data taken after this change-over. For the 
most part, this step in frequency has little effect on the 
data as a whole. Nevertheless, fields observed during the 
frequency change were reduced and stacked separately to 
test for any possible systematic effects introduced by the 
frequency change. No significant discrepancies were found 
(section [3.51 . 

Finally, changes in system temperature and in atmo- 
spheric opacity with elevation are corrected, based on the 
monitoring of modulated noise signals injected into the VSA 
system. Each antenna has its own noise diode to monitor 
changes in gain or Tsys for each antenna. This is particu- 
larly useful for identifying periods of bad weather and other 
problems with individual antennae. The correction is typi- 
cally a few per cent. Observations with large correction fac- 
tors (> 20 per cent), usually caused by bad weather or warm 
receivers, are discarded. 

3.4 Absolute flux calibration 

The absolute flux calibration of VSA data is determined 
from observations of Jupiter assuming a brightness temper- 
ature, Tjup. Earlier VSA results assumed r,jup = (152±5) K 
(3 per cent accuracy in temperature) at 32 GHz, as reported 
in Mason et al. (1999). The recent 1-year WMAP data (Ben- 
nett et al. 2003a) are scaled to their own measurement of the 
CMB dipole with an accuracy of 0.5 per cent (Hinshaw et 
al. 2003b). Page et al. (2003) have measured Jupiter to be 
Tjup = (146.6 ± 2.0) K at 33.0 GHz corresponding to an 
accuracy of 1.5 per cent in temperature terms, or equiva- 
lently 3 per cent in the CMB power spectrum (AT^). Note 
that the WMAP error for T!jup is dominated by beam shape 
uncertainties in the WMAP 1-year analysis. 

We have used the new WMAP temperature for Jupiter 
as the basis for our absolute temperature scale (Table 
In principle, this reduces the VSA calibration error in the 
power spectrum from 7 to 3 per cent. Indeed, the new value 
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off-diagonal terms remain, they will tend to average away 
to zero due their complex (real and imaginary) nature, as 
more days of data are combined. A further test was to com- 
pare the r.m.s. level well outside the primary beam with 
that computed for an "auto-subtracted" map (which is free 
from CMB or any other correlated signal; see Paper II). We 
found that the noise is indeed uncorrelated from visibility to 
visibility which assures us that the noise really is decreasing 
with the square root of time. 

The reduced data are stacked and combined in a va- 
riety of ways to look for possible contamination or other 
instrumental affects. By stacking the complete data set in 
Hour Angle, for a given field, any spurious signals (see Pa- 
per II) that had not been filtered out, would add coherently 
for each baseline and were therefore easy to detect and flag. 
For the extended array, very few instances of this effect were 
seen. In addition, the real and imaginary components of the 
visibility stacks were inspected in the u, v plane. 

We note that where possible, the data are reduced in 
parallel, by at least two of the three institutions in the VSA 
collaboration, in order to perform consistency checking. We 
compared the fields analysed independently by different in- 
stitutions to check that there were no major data reduction 
issues due to variations in cut-off levels, filtering parameters 
and subjective flagging of the data. No significant differences 
were found. The same comparisons were made in the power 
spectra which were also found to be completely consistent 
with each other. 

As an additional check on data consistency, the data for 
each VSA field were split by epoch and the statistic was 
computed. Since the sky signal is the same in both halves 
of the data, the statistic can be used to test for the pres- 
ence of systematic errors. Table 2] gives the values and 
associated significances. The term significance in this con- 
text represents the probability of exceeding the observed 
value in the cumulative distribution function. For each field, 
and averaged over all fields, the reduced is ~ 1 with a 
combined significance of 0.50. This indicates that the splits 
are statistically consistent with each other. A wide range of 
non-Gaussianity tests were performed on individual obser- 
vations and on the stacks. The techniques and results from 
earlier analyses are presented in Savage et al. (2004) and 
Smith et al. (2004). Detections of non-Gaussianity in the fi- 
nal stacked data were traced to a small portion of the data 
and hence cannot be cosmological in origin. These data were 
excised before calculating the final power spectrum in sec- 
tion |^ 



4 FOREGROUNDS 
4.1 Source subtraction 

4-1.1 Observational strategy 

Discrete sources are the largest foreground for the VSA. The 
power spectrum of Poisson-distributed sources increases as 
f (Taylor et al. 2001) while that of the CMB decreases ex- 
ponentially with increasing I. The point-source contribution 
will therefore dominate at higher ^-values in the absence of 
an effective source-subtraction strategy. We carefully choose 
regions of the sky to not contain bright (^ 500 mjy at 
1.4 GHz) radio sources f section 13. H when extrapolated to 



Table 4. The values for data splits on each of the VSA fields. 
In each case the visibility data from each field were split in two 
according to epoch and the of the difference vectors formed. 
Also tabulated are the number of degrees of freedom (d.o.f.) and 
the significance of each value (see text). Fields which strad- 
dle the observing frequency change-over date are separated into 
34 GHz ('A') and 33 GHz ('B') stacks. 



Field 


d.o.f. 




Significance 


VSAlE 


4778 


4851 


0.23 


VSAlF 


3871 


3964 


0.14 


VSAIG 


4070 


3752 


0.99 


■\ T'C A 1 XJ A 

VbAlH_A 


3867 


3790 


0.81 


T A 1 TJ TD 


3655 


3800 


0.04 


VSAIJ 


3845 


3890 


0.30 


■\ T'Q A 1 Ti^ 

V oAilv 


3706 


3691 


0.56 


VSAlL 


3742 


3733 


0.54 


■\ TC A OTP 


5002 


5170 


0.05 




3831 


3763 


0.78 




4314 


4297 


0.58 




4142 


4052 


0.84 


■\ TC A O T A 


4120 


4161 


0.32 


■\ 7C A O T TZ> 

VbA2J_r> 


3963 


4142 


0.02 


VbAiK 


3827 


3731 


0.86 


VSA2L 


3820 


3996 


0.02 


■\ re A OTT' 


4937 


5011 


0.22 


\ Td A OT7 

VbAdr 


4300 


4325 


0.39 


■\ 7c A or^ 
VbAdCj 


4970 


5060 


0.18 


■\ T'C A OXJ A 

VbAoH_A 


4236 


4252 


0.42 


VSA3H_B 


3926 


3892 


0.64 


VSA3J 


4105 


4081 


0.60 




OfUO 




U.Ol 


VSA3L 


4354 


4338 


0.56 


VSA5E 


3692 


3427 


0.99 


VSA5F 


1627 


1651 


0.33 


VSA5G 


3271 


4762 


0.00 


VSA6E.A 


701 


669 


0.80 


VSA6E.B 


1297 


1279 


0.63 


VSA6F 


3578 


3507 


0.80 


VSA6G 


2085 


2119 


0.29 


VSA7E.A 


711 


716 


0.43 


VSA7E.B 


2623 


2572 


0.75 


VSA7F 


3560 


3532 


0.63 


VSA7G 


3625 


3639 


0.43 


VSA8E.A 


3039 


3015 


0.62 


VSA8E.B 


1827 


1796 


0.69 


VSA8F 


2003 


1986 


0.60 


VSA8G 


3983 


4064 


0.18 






Mean = 


0.50 



33 GHz from NVSS (1.4 GHz) and GB6 (4.85 GHz) cat- 
alogues. However, the contamination from sources below 
~ 500 mJy is still significant, particularly at the highest 
^-values. 

To correct for any sources unidentified in the above sur- 
veys, the VSA employs a unique two-stage source subtrac- 
tion strategy. The first stage is to survey the regions ob- 
served by the VSA with the Ryle Telescope at 15 GHz to 
a limiting flux density of ~ 10 mJy (Waldram et al. 2003). 
This limit was chosen primarily to make sure that all sources 
that could potentially contaminate VSA data are located, 
even if they have rising spectra, a < 0, between 15 and 
30 GHz. Only the very steepest inverted spectrum {a < — 1) 
sources will be missed using this technique. The results from 
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the recent 9C survey (fig. 9 in Waldram et al. 2003), from 
the Ryle Telescope at 15 GHz and a completeness limit of 
25 mjy, indicate that less than 1 per cent of their radio 
sources have a < —1. We therefore expect little or no effect 
from such a rare population of sources at these flux density 
levels. The second stage, is to follow-up these sources with a 
single-baseline interferometer at the same frequency as the 
VSA and at the same epoch; each source is observed many 
times with an interval of a few days. The source subtrac- 
tion interferometer consists of two 3.7 m dishes with a 9 m 
N-S separation corresponding to a resolution of ~ 3 arcmin 
and has a sensitivity of « 340 mJy s~^/'^. The dishes are 
located in separate identical enclosures to minimise ground 
spill-over and cross-talk. 

4-1-2 Source subtractor observations and data reduction 

For a typical observation, the drive system of the source 
subtractor takes the coordinates of the previously identified 
15 GHz sources within a radius of ~ 2° of the VSA field 
centre. A different subset (typically ~ 30) of these sources is 
monitored over the course of each ~ 5-hour observing run. 
In this manner, the fiux density of the sources is sampled 
over the course of the observations of a given field. This takes 
into account source variability on time-scales longer than a 
few days. 

The single baseline interferometer has identical back- 
end hardware to those of the main array and the data are 
processed in a similar way. The primary fiux calibrator for 
the source subtractor is the planetary nebula NGC 7027. We 
assume a fiux density of (5.45 ± 0.20) Jy at 32.0 GHz and 
a spectral index a = (0.1 ± 0.1) (Mason et al. 1999), cor- 
responding to a 4 per cent calibration uncertainty. Vari- 
abUity of NGC 7027 is below 3 per cent (Peng et al. 2000) 
at VSA frequencies, apart from a secular decrease with time 
of ~ 0.6 per cent per year (Ott et al. 1994). Phase calibra- 
tion is provided by interleaved observations of the brightest 
sources {jt 100 mJy) in each field. 

4-1-3 Results 

The result of the source monitoring programme is two-fold. 
Firstly, the 33 GHz flux densities of the sources in the VSA 
flelds are known (at the time of the main array observation) 
and are subtracted from the visibility data down to a level 
of 20 mJy. The details of the complete source survey will be 
presented by Cleary et al. (in prep.). Secondly, the source 
counts derived from the observations can be used to apply 
a statistical correction due to fainter sources below the sub- 
traction limit. This analysis assumes that the contribution 
from clustered sources is negligible compared to that from 
the Poisson component. This is a reasonable assumption on 
the basis of estimates for the Planck 30 GHz channel by 
Toffolatti et al. (1998). 

Source monitoring for VSA source subtraction was per- 
formed at 34 GHz and 33 GHz. In order to derive the source 
counts, the flux density of each source measured at 34 GHz 
was corrected to 33 GHz using its af5.2 spectral index. An 
average was then taken of the source measurements over 
the entire period of extended array observations. Of the 
453 sources monitored, 131 were found to be > 20 mJy at 
33 GHz. 



The 33 GHz source counts at flux densities less than 
~ 5 mJy are unknown; however the counts are expected to 
flatten with decreasing flux density. Hence, simply extrap- 
olating the measured VSA source counts to fainter fluxes 
would over-estimate the residual source contribution. In the 
absence of low fiux density data, the measured VSA source 
counts were used to re-scale the Toffolatti et al. (1998) 
33 GHz differential source count model. Integrating the 
rescaled model up to 20 mJy, we get an estimate of the 
contribution from faint sources below the source subtrac- 
tion limit. The residual source power spectrum (~ 210 /iK^ 
at £ = 1000) is then binned using the VSA window functions 
and directly subtracted from the band-power estimates as an 
uncorrelated statistical correction. The correction applied to 
these data is AT^^^ = 210 x (^/lOOO)^ in units of fiK^. 



4.2 Galactic foregrounds 

4-2-1 Overview 

At frequencies of 30 GHz, particularly at large angular 
scales {£ < 100), emission from the Galaxy can contami- 
nate CMB data (Banday et al. 2003; Bennett et al. 2003b). 
There are currently three well-established diffuse Galactic 
foregrounds: synchrotron emission from relativistic electrons 
spiralling in the Galactic magnetic field, free-free (thermal 
bremsstrahlung) emission from ionized gas, and vibrational 
dust emission. The power spectrum of Galactic emission falls 
with increasing £ (for example, see Giardino et al. 2001) 
and hence it is not expected to be a major contaminant at 
smaller angular scales. The VSA is insensitive to large angu- 
lar scales (^ 1°) due to the incomplete u, v coverage and pri- 
mary beam attenuation. Furthermore, the VSA fields have 
been carefully selected to be located in "cold" areas of sky 
f sect ion |3. in . at high Galactic latitudes, thus minimising any 
potential contamination from the Galaxy. There is however, 
evidence for another foreground component, recently nick- 
named "foreground X" (de Oliveira-Costa et al. 2004), that 
has been shown to be strongly correlated with FIR emis- 
sion (A ~ 100 i-im) in the frequency range ~ 10 — 50 GHz 
(Kogut et al. 1996; Leitch et al. 1997; de Oliveira-Costa et 
al. 2002; Mukherjee et al. 2002; Banday et al. 2003; Bennett 
et al. 2003b). A recent discussion on the possible origins of 
the anomalous emission, including spinning dust (Draine & 
Lazarian, 1998) and magnetic dust grains (Draine & Lazar- 
ian, 1999) is given by Banday et al. (2003). The 1-year 
WMAP data also detect the dust-correlated component in 
the K, Ka and Q bands (Bennett et al. 2003b). Using the 
statistical distribution of derived spectral indices, they in- 
terpret this to be synchrotron emission with a relatively fiat 
spectral index, /3 ~ 2.5 (T oc «^~'^), and constrain the spin- 
ning dust component to contribute less than 5 per cent to 
the foregrounds. They argue that the 408 MHz map is not 
a good template for synchrotron emission at higher frequen- 
cies (J; 20 GHz) because it is dominated by steep spectrum 
synchrotron emission with (3 ~ 2.9. This is a contentious 
issue which has already been re-examined but with no con- 
sensus on the origin of the anomalous component (e.g. La- 
gache 2003; Banday et al. 2003; Finkbeiner et al. 2002; 
Finkbeiner 2003a; Casassus et al. 2004; de Oliveira-Costa 
et al. 2004). There is nonetheless clear evidence for a dust- 
correlated component, at frequencies of ~ 30 GHz, with a 



8 C. Dickinson et al. 



typical coupling coefficient to the SFD98 100 /xm map of 
Tb//ioo ~ 10 ^K/(MJy sr"^). Indeed, the VSA has observed 
selected dust clouds to further investigate the Galactic emis- 
sions at 33 GHz (Dickinson et al., in prep.). 

4.2.2 Galactic foreground estimates 

To quantify the Galactic foreground we use external tem- 
plates since a spectral analysis is not possible for single fre- 
quency data. The template maps are (i) the 408 MHz all-sky 
map at a resolution of 51 arcmin (Haslam et al. 1981) for 
steep spectrum synchrotron, (ii) Ha data from the Wiscon- 
sin H-Alpha Mapper (WHAM) at 1° resolution (Haffner et 
al. 2003) for free-free emission and (iii) the 100 ^m map at 
a resolution of 6.1 arcmin (SFD98) for dust-correlated emis- 
sion. The synchrotron and free-free templates do not have 
adequate resolution to permit a full cross-correlation anal- 
ysis, although high resolution Ha data will be available in 
the near future (Dennison, priv. comm.) from the Virginia- 
Tech Spectral-line Survey (VTSS; Dennison, Simonetti & 
Topasna 1998). We therefore calculate the r.m.s. fluctua- 
tions in each region at the resolution of 1° in a similar man- 
ner to the analysis used in Paper H using simple foreground 
models to convert to temperature fluctuations at 33 GHz. 
We use circular regions over a 4°.5 and 3°.0 diameter for the 
7-field mosaics and 3-field mosaics respectively. We assume 
a spectral index for synchrotron oi 13 = 2.9 and a conver- 
sion factor from Ha intensity, in units of Rayleigh (where 
IR = 10®/47r photons s~^ m~^ sr^^), to brightness tem- 
perature at 33.0 GHz, of 5.0 R^^ (Dickinson, Davies & 
Davis 2003; Finkbeiner 2003b). The r.m.s. values were neg- 
ligible (;$ 10 ^jK^) in the VSA fields compared to the CMB 
fluctuations (> 1000 ^JK^). 

For the dust-correlated component, for which we have 
a high resolution (6.1 arcmin) template, we smoothed the 
SFD98 100 /^m map to 22 arcmin {(. ~ 1000) and estimated 
the foreground levels by assuming a typical coupling coef- 
flcient between brightness temperature at 33 GHz and the 
100 /xm intensity of Tb//ioo = 10 /iK/(MJy sr"^). The r.m.s. 
estimates are listed in Table |H1 At ^ = 1000, the CMB fluc- 
tuations are at ATj^n^^ ~ 1000 /iK^ while the r.m.s. fore- 
ground estimations are H 10 ^K^. Hence, for most of the 
VSA regions, the Galactic emission is essentially negligi- 
ble unless they have an unusually large coupling coefficient 
(^ 30 /iK/(MJy sr~^)). From the r.m.s. estimates in Table |3 
the VSAl and VSA5 regions seem to have stronger Galactic 
emission compared to the other flelds. Nevertheless, these 
are still below the CMB anisotropies. Furthermore, there 
is no evidence of peculiarities when comparing these flelds 
with the other VSA flelds both in terms of the reconstructed 
maps f section l5.1^ or the power spectra of individual VSA 
regions discussed further in section]^ 

We also estimate the foreground power spectra in the 
VSA regions by simulating observations of the 100 /xm map 
in the VSA flelds based on the actual u, v positions of the 
visibilities and taking into account the primary beam. The 
resulting power spectra are well-modelled by a simple power- 
law function in each case and the expected fall-off in power 
with increasing £ agreed with previous analyses (SFD98), 
with power-law indices in the range —0.5 to —1.0. Of course, 
the exact scaling law depends on frequency and position on 
the sky. We fltted VSA data to the 100 fim template for each 



Table 5. Galactic dust-correlated foreground power estimates 
(r.m.s) at £ = 1000 for each VSA region based on the SFD98 
100 /im map. We assume a correlation coefficient, between the 
foreground map (in units of MJy sr~^) and 33 GHz brightness 
temperature (in units of /xK), of 10 /xK/(MJy sr"-*^). 



VSA region 


r.m.s. power fluctuations 




(/.K2) 


VSA 1 


31.9 


VSA 2 


1.3 


VSA 3 


7.9 


VSA 5 


88.4 


VSA 6 


9.6 


VSA 7 


2.3 


VSA 8 


9.0 



individual VSA fleld. This was difficult due to low-signal-to 
noise ratios and the relatively stronger CMB anisotropies. 
The fltted values were in the range — 30 fiK/{MJy st~^) 
but with large errors and no correction for chance align- 
ments of the template maps and the CMB. Taking the av- 
erage over all flelds was very close to the canonical value of 
10 MK/(MJy sr"^). Since we are only interested in the over- 
all effect in the combined power spectrum, we chose to stick 
with the canonical value. It was used to re-scale the power 
spectra fits to each region and hence provide a statistical 
correction to the final VSA power spectrum in section IF!^ 
The scatter on the derived coupling coefficients both here 
and from other experiments, at similar frequencies and an- 
gular scales, suggests that the overall scaling could be wrong 
by up to a factor of ~ 2 at most. However, the correction is 
very small for these data and hence such an error would have 
only a minor effect at the lowest ^-values {£ 100) where 
the VSA has little or no sensitivity; this is discussed further 
in section |H] The statistical correction applied to the final 
power spectrum has the form ATq^i = 4.1 x 10^ x £~°-''^. 
This corresponds to 28 /iK^ at £ = 1000. A more detailed 
analysis of Galactic foregrounds will be described in a forth- 
coming paper (Dickinson et al., in prep.). 

4.3 Galaxy clusters foreground — the SZ effect 

Galaxy clusters can be a serious contaminant for CMB ob- 
servations. Large reservoirs of hot gas (T ~ 10^ K) in the 
intracluster medium re-scatter the CMB photons by inverse 
Compton scattering. Since the optical depth is typically 
~ 1 per cent, this creates a small distortion of the CMB 
frequency spectrum — the Sunyaev-Zeldovich Effect (Zel- 
dovich & Sunyaev 1969; Sunyaev & Zeldovich 1970), hence- 
forth SZE. For frequencies below the CMB spectrum peak at 
217 GHz, the SZE corresponds to a temperature decrement 
along the line of sight to the cluster. The largest clusters, at 
relatively low redshift {z < 0.1), can easily be detected by 
the VSA (Rusholme 2001; Lancaster et al. 2004) and hence 
it is crucial to avoid these when making CMB observations. 
Fortunately, the majority of nearby clusters have been iden- 
tified by optical surveys such as the Abell catalogue (Abell, 
1958). The VSA regions are located away from known clus- 
ters fsecti0n r3.ll . 

One potential contaminant, which we have not dis- 
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cussed so far, is the Poisson distributed contribution due 
to SZE decrements from virialized clusters at higher red- 
shifts. For 2000 < I < 4000 this is likely to have the same 
^-dependence as the point source contribution, albeit with 
a lower amplitude. If normalisation cts = 0.9 then we esti- 
mate, using the cluster model described in Battye & Weller 
(2003), a confusion noise of ~ 1 mjy beam~^, well below 
the noise level in typical VSA maps, of ~ 6 mJy beam~^. If 
the excess power detected by CBI (Mason et al. 2003; Read- 
head et al. 2004) is due to the SZE then ag ~ 1.1 (Komatsu 
& Seljak 2002) and we estimate 2.5 mJy beam^^ compati- 
ble with an extrapolation of the observed power spectrum 
from I = 3000. In either case, the SZE foreground is not 
significant for current VSA data. 



5 RESULTS 

5.1 Maps 

The primary goal of the VSA is to measure the CMB power 
spectrum over a wide range of angular scales. Sky maps also 
provide an important and complementary view of the data. 
To create maps of the CMB we used a maximum entropy 
method (MEM) described in Maisinger, Hobson & Lasenby 
(1997). The input data are the same binned visibility data 
used for deriving the power spectrum in section 15.21 We 
chose to use the "memuv" option which reconstructs the 
Fourier modes in the u, w-plane; these are then Fourier trans- 
formed to make maps. To test the robustness of the map- 
making software, we tried a range of options with various 
parameters when producing the maps. We found no signif- 
icant differences between these maps. The MEM algorithm 
is preferred over other deconvolution algorithms since it is 
optimized for extended emission by assuming a flat sky as a 
prior. The CLEAN algorithm (Hogbom, 1974), on the other 
hand, is optimized for unresolved sources which are treated 
as delta functions convolved with the synthesized beam. 

The MEM-based CMB maps made from VSA data in 
the extended array, are presented in Fig. 13 Foreground radio 
sources have been removed as explained in section llTTI CMB 
anisotropies are clearly detected in all 7 VSA regions. 

The primary beam envelope of an interferometer per- 
mits a separation of the astronomical signal and the noise. 
Within the primary beam, the map contains both CMB 
and noise components while outside the primary beam, only 
noise remains. This is not the case for the MEM maps in 
Fig. 13 so we also created maps based on the CLEAN de- 
convolution algorithm with the IMAGR task in the AIPS 
package. Table|S|lists the estimates of CMB signal and ther- 
mal noise for individual VSA fields. Significant variations in 
the signal-to- noise ratio (SNR) are evident (SNR ~ 1 — 3), 
primarily due to the varying integration times after flagging 
and filtering of the data. Note that the signal-to-noise ratio 
in the mosaiced maps (Fig. |3 is better due to the overlap- 
ping of the individual fields. 

5.2 Power spectrum 

The final visibility data are binned into square w, v cells, each 
9 wavelengths on a side, to oversample the data while reduc- 
ing the number of data points by a factor of ^ 1000. Sources 



Table 6. Comparison of r.m.s. of the signal and thermal noise, 
all units mJy beam~^, for each VSA field, measured from the 
CLEANed maps. The r.m.s. noise level in the centre of each map 
and the corresponding residual CMB signal observed in each field 
are given. Fields with 2 values (A and B) are for data at 2 fre- 
quencies (34 and 33 GHz respectively). 
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are subtracted using position and flux density information 
as described in section im This results in 33 visibility files, 
one for each VSA pointing, each with ~ 10^ data points. 

The binned visibilities form the basic input to the max- 
imum likelihood analysis for the CMB power spectrum. 
We used the Microwave Anisotropy Dataset Computational 
sOftWare (MADCOW) described by Hobson & Maisinger 
(2002). 

The main binning of the VSA data reduction was cho- 
sen to Nyquist sample and to approximately match the ex- 
pected phase of the peaks and troughs in the power spec- 
trum as closely as possible. This means that there is a bin 
corresponding to the ^-value expected for each peak in the 
power spectrum and an additional bin between the peaks. 
Note that this approach does not bias the parameter estima- 
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Figure 4. VSA window functions plotted as a function of £ for 
the main binning. Each curve is a single bin out of a total of 16 
bins. Curves are solid and dashed line for clarity. The window 
functions are, in general, well-modelled by Gaussian functions. 



array data sets, are listed in Table |H| for two independent 
binnings'^. The final combined VSA power spectrum is pre- 
sented in Fig.|S] Bright foreground sources (section k.H have 
been subtracted from the visibility data before power spec- 
trum estimation. The small corrections for residual sources 
(section 14. 1.311 and diffuse foregrounds (section 14.2.21 have 
been applied and are also plotted in Fig. |^ The error bars 
were calculated from the probability likelihood functions by 
enclosing 68 per cent of the area centred on £h, the weighted 
median £ value for each bin. Calibration uncertainty (3 per 
cent; section lX^ is not included. Sample variance is included 
in the error estimates. Note that the extended array data 
have little sensitivity at ^ J; 300 and hence the first 3 bins 
come mainly from data taken in the compact array described 
in Paper III. 

The VSA data are also plotted in Fig. |S| along with 
recent data from other CMB experiments; WMAP (Bennett 
et al. 2003a), ACBAR (Kuo et al. 2004) and CBI (Readhead 
et al. 2004). Only a single binning of the data shown. No 
corrections are made to the external measurements for the 
binning schemes; they are plotted directly. Nevertheless, the 
overall consistency between these experiments is remarkable. 



Table 7. The correlation matrix for the combined VSA data 
(main binning only), normalised to 1.0. The values of the matrix 
for which dj is not reported, can be assumed to be zero. The 
final column gives the diagonal elements of the covariance matrix 
in units of lO'' X /^K*. 
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tion towards the concordance model (see for example, Kuo 
et al. 2004) — it merely acts as a matched filter that tries 
to extract maximum information from the fluctuations in 
the power spectrum. The bin positions and widths are the 
same as those used in Grainge et al. (2003). The consider- 
ably larger contiguous sky area in the seven fleld mosaics 
results in reduced correlation between adjacent bins. The 
inverse correlation matrix for the combined VSA data set 
is given in Table |7| The window functions are calculated as 
described in Paper III and are plotted in Fig. 2] they are 
well-fltted by Gaussian functions. In principle, finer binning 
is possible but at the cost of increasing the correlations and 
reducing the signal-to-noise ratio. 

The flnal band powers calculated from the complete 
VSA data set, including both the compact and extended 



5.3 Morphological power spectrum analysis 

The shape of the power spectrum was analysed following 
the technique described in Odman et al. (2003). By fitting 
multiple Gaussians, it calculates the best-fit to the data. 
Each Gaussian is defined by its amplitude in temperature 
AT, centre multipole I and width a. 

The data were split into subsets according to the nature 
of each experiment. The low-frequency (LF) subset consists 
of DASI (Halverson et al. 2002), CBI (Mason et al. 2003; 
Pearson et al. 2003) and VSA data (this work) , with frequen- 
cies ~ 30 GHz. The high frequency (HF) subset, correspond- 
ing to frequencies of ^ 100 GHz, comprises Archeops (Benoit 
et al. 2003), ACBAR (Kuo et al. 2004), BOOMERanG (Net- 
terfield et al. 2002) and MAXIMA (Lee et al. 2001) data. 
An additional analysis on the VSA data alone was carried 
out to check for consistency. WMAP data were not included 
to simplify the low and high frequency split. 

The results of this analysis — the amplitudes and centre 
^- values for each peak in the spectrum — are in Table |5] LF 
data do not tightly constrain the first peak position owing 
to the lack of data points at low I. This appears as a double 
peak in the likelihood for l\ and ATi . The uncertainty on the 
first peak does not propagate at higher multipoles where the 
secondary peaks are increasingly well-defined. HF data pro- 
vide constraints on the 1st, 2nd and 3rd peaks whereas LF 
data provide constraints on the 2nd to the 5th peaks. An im- 
portant aspect of this analysis is the apparent discrepancy, 
at ~ 2(7 level, between LF and HF data on the amplitude of 
the 3rd peak. The LF data prefer a slightly higher 3rd peak 
of ~ 48 /iK, compared to the HF data which prefer a lower 
value of « 40 /iK. This feature has been noted before and 
mostly accounts for the high baryon fraction derived from 
LF data alone. 

^ The VSA data are available at the following URL: 
|http : //wHW. jb.man. ac ■uk/research/vsa/vsa_results .html| 
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Figure 5. The CMB power spectrum as measured by the VSA by combining the extended and compact array data for all 7 VSA regions. 
The errors represent Icr limits from both thermal noise and sample variance. Two alternate binnings (filled red circles and unfilled blue 
circles) are displayed fTablelSl. Absolute calibration is accurate to 3 per cent and is not included in the errors. The bright radio sources 
have been subtracted. The smaller corrections for residual sources and Galactic foregrounds have been made and are shown as a dashed 
line and dotted line respectively. 
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Figure 6. A comparison of the VSA data (red filled circles) with those from WMAP 1-year (crosses), CBI (filled blue squares) and 
ACBAR (filled stars). All errors are at the Icr level. There is good agreement between the different experiments. 



Table 9. Results of fitting Gaussians to the CMB power spectra as measured by various CMB experiments. Here we show four combi- 
nations of data (see text), including low frequency (LP; !^ ~ 30 GHz) and high frequency (HP; u ~ 100 GHz) data. Amplitudes are in 
units of /iK. Errors are at the 68 per cent confidence level; "nc" is where there are no constraints. 
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Table 8. The CMB band-powers (in /^K^) for main and offset 
binnings of the complete VSA data set combining both compact 
and extended array data. The quoted ^-ranges are nominal bin 
limits when the window functions fall to half the peak value, while 
Ih is the weighted median value of the relevant window function. 
The error bars enclose 68 percent of the likelihood and include 
sample variance. 
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This type of model-dependent analysis is often difficult 
to interpret. Specifically, the 4th and 5th peaks are not real 
detections and are merely a consequence of how the fitting 
procedure accounts for the exponential drop in power due 
to the damping tail. Nonetheless, the structure of the CMB 
power spectrum, at least at ^ < 1000, is relatively well- 
defined. There is good agreement between this analysis and 
a similar one made for the 1-year WMAP data on the 1st 
and 2nd peaks (Page et al. 2003). 



6 DISCUSSION 
6.1 Maps 



The largest foreground component for the VSA is dis- 
crete sources. The stronger sources can be identified in the 
raw VSA maps, before source subtraction. They show up 
as large positive temperatures in the maps, although they 
can coincide with negative CMB features (or noise), mak- 
ing faint sources difficult to identify in the sky maps. A 
comparison of an un-subtracted map with one that has had 
sources removed is presented in Fig. |7| These maps have 
been deconvolved with the CLEAN algorithm to optimise 
detection of discrete sources. In VSA6E, there are 3 rela- 
tively strong sources, even after correction for the primary 
beam attenuation, with apparent fiux densities of 91, 84 and 
39 mjy. The r.m.s. noise in the data is ~ 10 mjy. The re- 
quirement for accurate source subtraction, as used by the 
VSA, is evident. The VSA maps allow a comparison to be 
made with other CMB data. The recent WMAP data re- 
lease (Bennett et al. 2003) has provided 5 all-sky maps at 
frequencies centred at 22.8 (K-band), 33 (Ka-band), 40.7 
(Q-band), 60.8 (V-band) and 93.5 GHz (W-band) with res- 
olutions ranging from 49.2 arcmin (K-band) to 12.6 arcmin 
(W-band). The signal-to- noise ratio of WMAP data at the 
VSA resolution is ~ 1 and hence much of the CMB signal 
is lost in the noise. The actual noise level in the WMAP 
data depends on position due to the scanning strategy of 
the WMAP sateUite. For the 1-year WMAP data release, 
the noise is ~ 100 — 200 /iK per 12.6 arcmin pixel in the 
VSA regions, compared to ~ 25 /iK in the VSA mosaiced 
maps. 

The WMAP team provide their own foreground-cleaned 
CMB maps by combining all the WMAP channels and opti- 
mising the individual weights to a CMB signal. However, the 
resolution is constrained by the K-band channel at a resolu- 
tion of 49 arcmin. More recently, Tegmark, de Oliveira-Costa 
& Hamilton (2003) have produced a foreground-cleaned 
CMB map at the resolution of the highest frequency channel 
(12.6 arcmin). In Fig.|H] we compare the VSA data for region 
VSAl, with the foreground-cleaned CMB map of Tegmark 
et al. (2003), based on WMAP 1-year data at a resolution 
of 12.6 arcmin. The colour scale uses the same range ( — 200 
to -1-200 /iK) as that used in Fig. |3 Some features, both 
positive and negative, can be seen on both maps. The cor- 
relation is however restricted due to the larger noise level 
in the WMAP map for this field. The minimum and maxi- 
mum values for these maps are typically ~ ±400 ^K, while 
for the VSA maps the minimum and maximum values are 
~ 250 /xK, in agreement with simulated CMB maps based 
on the concordance model and a noise level of ~ 100 /xK. 
Fig. |51 also shows the SFD98 100 /^m dust map, smoothed 
to 11 arcmin resolution. For the VSAl region, which con- 
tains moderate dust emission, there are no obvious Galactic 
signals in either the VSA or WMAP data for this region. 
It is evident that the strongest Galactic emissions are on 
large angular scales {i < 100), at least at high Galactic lat- 
itudes (l&l ^ 30°). In general, the VSA maps are consistent 
with each other with no obvious sign of contamination from 
foregrounds or other correlated signals. 



An important application of CMB maps is to look for the 
signature of any non-Gaussian structures which may be fore- 
grounds or primordial in origin, such as cosmic strings which 
would appear as a discontinuity on a map, or a SZE which 
would produce a negative temperature against the CMB. 



6.2 Power spectrum 

The VSA power spectrum (Fig. |SJ clearly shows the ex- 
istence of the first three acoustic peaks and the exponen- 
tial fall-off in power towards higher I. The WMAP all-sky 
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Figure 7. VSA maps of VSA6E field before and after source subtraction (Cleary 2003). The maps have been deconvolved using the 
CLEAN algorithm. The FWHM of the primary beam is indicated by the circle. Astronomical signals are inside the primary beam while 
the edges of the map, outside the primary beam, contain noise. The crosses indicate the positions of the three brightest sources in this 
field (see text). The r.m.s. noise level is 10 mjy. 




Figure 8. Comparison of maps in region VSAl. Left: Foreground-cleaned CMB map based on WMAP 1-year data as re-analysed 
by Tegmark et al. (2003), with a resolution of 12.6 arcmin. Middle: Source-subtracted CMB map from VSA data at a resolution of 
~ 11 arcmin. Right: SFD98 100 ;jm map (in units of MJy sr~l) smoothed to 11 arcmin. The WMAP noise level is ~ 100 /iK (lo"). 



data provides measurements limited by cosmic variance, at 
I < 350. At moderate to high £ values {I 800) there are 
fewer measurements of the CMB power spectrum. Recently, 
the ACBAR and CBI instruments have released new results 
which extend the ^-range further. The morphological anal- 
ysis presented in section 1?^. 81 also suggests very good agree- 
ment between experiments, particularly on the 1st and 2nd 
peaks. The slight discrepancy in the amplitude of the third 
peak might suggest a change in the cosmological model. 
However, since the discrepancy is between low and high fre- 
quency data, it is possible that this could be a consequence 
of some systematic effects either in the LF or HF data. 

An important aspect of these results is the improved 
absolute calibration. We have used the WMAP brightness 
temperature for Jupiter resulting in a reduction of the VSA 
values by 8 per cent and an absolute uncertainty of 3 per cent 



in the power spectrum. Interestingly, with the new calibra- 
tion scheme, the VSA data are slightly lower in amplitude, 
at least over the ^-range 350 — 700. In fact, if we had used 
the original calibration scheme, we would be in closer agree- 
ment to the 1-year WMAP power spectrum. We find that 
the best-fitting scaling factor is ~ 1.05 i.e. the new VSA 
data are lower than the WMAP data by 1.6a. 

Foregrounds in the VSA data are expected to have a 
negligible effect on the final power spectrum. The brighter 
sources have already been subtracted prior to power spec- 
trum estimation (section 14.111 . The statistical corrections for 
residual sources and the Galactic foreground are relatively 
small as depicted in Fig|3 Nevertheless, it is clear that the 
residual sources are becoming significant in the highest i- 
bins. This emphasises the need for deeper source surveys if 
we are to accurately measure CMB fiuctuations at £ 1500. 
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On the other hand, the diffuse Galactic foregrounds are 
stronger at low ^-values. In the ^-range covered here, the 
foreground maps in the VSA regions suggest typical Galac- 
tic signals of a few tens of /iK^ at £ — 1000 compared to the 
CMB at ~ 1000 fiK^. The exception to this is VSA5, and to 
a lesser extent VSAl. The location of these fields (see Fig.l^J 
explains why the foreground estimates (Table |^ are higher 
than in the other fields. The diffuse foregrounds may indeed 
be stronger (~ 100 fJ-K^) in these fields but they still are 
below the CMB anisotropies (~ 1000 /iK^). This is verified 
by the fact that there are no obvious foreground signals seen 
either in the VSA maps or in the individual power spectra 
for each field. Moreover, the relative weights of the different 
VSA fields means that the combined power spectrum would 
still remain virtually unaffected — VSA5 contributes just 5 
per cent to the total weight of the data as a whole. This 
is a major advantage of observing a number of independent 
fields before combining them to increase the signal-to-noise 
ratio in the power spectrum. 

As for other systematics, the data from each individ- 
ual field are tested using a Bayesian joint power spectrum 
and non-Gaussianity estimation discussed in Rocha et al. 
(2001). This method operates directly on the visibilities and, 
as such, is ideally suited to the detection of any residual 
systematics (Savage et al. 2004). Additionally, the method 
properly accounts for the presence of noise in the data, mak- 
ing it extremely sensitive to any non-Gaussian signal. How- 
ever, we find no evidence of significant residual systematics 
in the data. 



7 CONCLUSIONS 

We have presented high sensitivity and foreground sub- 
tracted measurements of the CMB power spectrum up to 
£ = 1500 observed with the VSA. The cosmological in- 
terpretation is described in a companion paper (Rebolo et 
al. 2004). The final extended array data set contains a factor 
of ~ 4 more than those presented by Grainge et al. (2003). 
The sky-coverage has increased by a factor of ~ 3 utilising 4 
additional 3-field mosaics and increasing the original 3-field 
mosaics to 7-field mosaics, resulting in a significant increase 
to the signal-to-noise ratio. Mosaicing techniques allow for 
smaller correlations between adjacent bins and hence, can 
give finer resolution in £. Ifowever, in this paper we have 
used the same binning scheme as Grainge et al. (2003) for 
comparison and to reduce the bin-bin correlations ; the £ 
resolution is A£ ~ 60. We adopt a new absolute calibration 
scheme based on Jupiter temperature measurements with 
WMAP to give a 3 per cent accurate absolute scale in the 
power spectrum. This reduces absolute values by 4 per cent 
in temperature or 8 per cent in the power spectrum com- 
pared to previous VSA data. 

The VSA power spectrum (Fig. is in good agree- 
ment with data from other CMB experiments. There are 
now good signal-to-noise detections of the first 3 peaks in 
the CMB power spectrum and a damping tail at high £ which 
continues to ^ ~ 1500. The different instruments have dif- 
ferent potential systematic problems and they cover a wide 
range of frequencies. The agreement therefore indicates that 
none of the experiments are seriously contaminated by fore- 
ground emission or other systematic effects (see also Grif- 



fiths & Lineweaver (2004)). From the analysis of the peak 
structure fsection l5.3ll . there is extremely good agreement on 
the peak structure between experiments, except for a slight 
discrepancy in the amplitude of the 3rd peak which may 
have important consequences for the cosmological models. 
The slight discrepancy (~ 1.6a) in the overall scaling be- 
tween VSA and WMAP may also be an important factor 
in terms of fitting cosmological models. This is discussed by 
Rebolo et al. (2004). 

In the near future, the VSA will be reconfigured with 
even larger horns and longer baselines. It will allow the VSA 
to increase the ^-range even further, up to a maximum of 
£ 3000, while maintaining good £ resolution without loss 
of overall temperature sensitivity. The assistance of high fre- 
quency radio source surveys, to identify and measure the fiux 
densities of sources over large areas of sky, will be a crucial 
role in keeping the discrete source foregrounds under con- 
trol. The One Centimetre Receiver Array (OCRA; Browne 
et al. 2000), a prototype multi-beam receiver system cur- 
rently begin tested on the Torun 32 m telescope, will be ideal 
for such purposes. Moreover, OCRA is currently operating 
at 30 GHz and hence direct source subtraction, without the 
need to extrapolate to other frequencies will be possible in 
the near future. 
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